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Abstract — The objective of this work is to develop a recursive, 
discrete time quantum filtering equation for a system that 
interacts with a probe, on which measurements are performed 
according to the Positive Operator Valued Measures (POVMs) 
framework. POVMs are the most general measurements one 
can make on a quantum system and although in principle 
they can be reformulated as projective measurements on larger 
spaces, for which filtering results exist, a direct treatment of 
POVMs can be more natural and simplify the filter computa- 
tions for some applications. Hence we formalize the notion of 
commuting (Davies) instruments which allows one to develop 
joint measurement statistics of two POVM type measurements 
as explicit functions of the POVMs. This allows us to prove 
the existence of a notion of conditional expectation POVM, 
which is essential for the development of a filtering equation. 
We demonstrate that under generally satisfied assumptions, the 
reduced model given by POVM elements is sufficient for the 
purpose of our quantum filtering task. 



I. Introduction 

The theory of filtering considers the estimation of the 
system state from noisy and/or partial observations (see 
e.g. [1]). For quantum systems, filtering theory was initiated 
in the 1980s by Belavkin in a series of papers [2], [3], 
[4], [5], Belavkin makes use of the operational formalism 
of Davies [6], which is a precursor to the theory of quantum 
filtering. He has also realized that due to the unavoidable 
back-action of quantum measurements, the theory of filtering 
plays a fundamental role in quantum feedback control (see 
e.g. [3], [5]). The theory of quantum filtering was indepen- 
dently developed in the physics community, particularly in 
the context of quantum optics, under the name of quantum 
trajectory theory [7], [8], [9], [10]. 

The basic model used to derive filtering equations for a 
quantum system uses a system-probe interaction. A quantum 
system, whose state needs to be estimated, is made to interact 
with a probe and the state of the system becomes entangled 
with that of the probe. After this interaction, an observable is 
measured on the probe and this measurement outcome is used 
to estimate the state of the system. The commutativity of 
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any system observable with any probe observable is used to 
develop a recursive Markov filtering equation for the system 
observables (see e.g. [11], [12] for an excellent tutorial 
review of these ideas). 

Suppose %s is the Hilbert space corresponding to the 
system whose state needs to be estimated and Hp is the 
Hilbert space of the probe. According to the classical von 
Neumann definition, any probe observable is a self-adjoint 
operator Q on Hp; the measurement of such an observable 
results in an outcome that is (stochastically) an eigenvalue 
of Q and the probe state after measurement gets projected 
onto the corresponding eigenspace of Q. As far as we are 
aware, all discussions on quantum filtering theory so far have 
assumed that the probe undergoes such a von Neumann mea- 
surement, also called projective measurement or Projection 
Valued Measure (PVM). However, a more modern treatment 
of quantum measurement theory shows that the most general 
possible quantum measurements that one can perform are 
the so-called Positive Operator Valued Measures 1 (POVMs), 
of which von Neumann measurements are merely a special 
case, where all the positive operators in the POVMs are 
commuting projections [6], [13]. See Section II for a brief 
overview of POVMs. 

POVMs on Hp can be reformulated as the restriction 
to Hp of a PVM on a larger space. However, there is 
no canonical PVM that corresponds uniquely to a given 
POVM. This is closely related to the fact that the state 
of a quantum system after a POVM measurement is not 
uniquely determined as a function of the POVM. To remedy 
the latter situation, Davies [6, Ch.3] has shown that one can 
associate a (non-unique) instrument to any POVM, which 
determines a completely positive map that specifies the 
state after measurement conditioned on the measurement 
outcome. However, there is again no canonical instrument 
that corresponds to a given POVM. Therefore, it is impos- 
sible to uniquely specify the post-measurement state for a 
given POVM measurement outcome unless the instrument 
associated with the measurement is known. As stated by 
Nielsen and Chuang [14, P. 91], "POVMs are best viewed 
as [...] providing the simplest means by which one can 
study general measurement statistics, without the necessity 
for knowing the post-measurement state." As such they are 
a minimal description of quantum measurements, so one 
can hope that the POVM formalism leads to more concise 
and fast filtering equations, suited for (possibly analog) 
implementation in real-time quantum feedback experiments. 

There are other reasons to develop a POVM-based filtering 

1 The terminology for such measurements is not standard and other terms 
such as Positive Operator Measures or generalized measurements are also 
used. 
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theory that shortcuts the lift to PVMs in larger spaces. For 
one, once the POVM theory is available it can be more natu- 
ral to use, as several practical measurement setups are based 
on POVMs, such as e.g. approximate position/momentum 
measurements [6, Ch. 3] or phase measurements [13], [15]. 
In some infinite-dimensional situations there can even be 
conceptual barriers to a PVM viewpoint. Indeed, with phase 
measurements after much investigation there is still no 
universal agreement on an acceptable PVM [16]. Finally, 
regarding system identification, the POVM associated with 
any experimental setup can (at least conceptually) be directly 
deduced from measurement outcomes; in contrast, in order 
to ascertain the associated instrument it is necessary to 
analyze the post-measurement state of the system (for more 
details see Section II-A). This is not generally feasible in 
practical experimental setups where often the measurements 
destroy the quantum state (e.g. photo-detection) and/or non- 
measurably alter it due to interaction with the environment. 

In this paper, we develop a discrete time filtering equation 
for the system state conditioned on POVM measurements 
performed on a probe. After reviewing the POVM formalism 
(Section II), we provide a general theory about the commut- 
ing instruments that are associated to POVMs (Section III). 
This allows us to define conditional expectations POVMs 
in the measurement-motivated spirit of [11], [12] (Section 
IV). In the setup consisting of a probe coupled to a target 
system, any (physically reasonable) instrument associated 
with a POVM acting only on the probe commutes with any 
instrument associated with a POVM acting only on the target 
system. In Section V we show how this allows to define a 
filtering equation for the system state conditioned on POVM 
measurement outcomes. A notable result is that this filtering 
equation is only a function of the probe POVM and does not 
depend on the associated instrument. 

II. Review: POVMS and associated instruments 

The POVM formalism is a standard part of most modern 
quantum information textbooks. We briefly review it in this 
section and refer the interested reader to [6, Ch. 3] for more 
details. 

Consider a quantum system with Hilbert space Ti, i.e. the 
system state is given by a density operator p which is 
a unit-trace nonnegative self-adjoint linear operator on T-L. 
We use * to denote the adjoint. Denote by C(T-L) the set 
of linear operators on Ti, by C+(H) C C(H) the set of 
self-adjoint nonnegative linear operators, and by S(Ti) C 
£+(H) the set of all possible density operators. Standard 
textbook treatment of quantum measurements assumes that 
any physically measurable quantity A is associated to a self- 
adjoint operator A : Ti — > Ti. Because A is self-adjoint, we 
have the spectral decomposition 2 



A = J2 



(1) 



where £1 is the set of eigenvalues of A and P u is the eigen- 
projection corresponding to eigenvalue u. Starting with a 

2 For clarity of explanation, we assume that A has a discrete spectrum. 
The discussion easily generalizes to the continuous spectrum situation. 



system in state p, according to von Neumann's measurement 
postulates we have: 

1) any measurement of the observable A gives some 
outcome lj G fi with probability Tr {pP w }, and 

2) after measurement outcome oj, the state of the system 
becomes 

9 ~ Tr {P uP Pu,}' 

The first postulate can be thought of as follows: the set 
fi is the set of all possible measurement outcomes of an 
experimental setup (A) and to each uj £ fl, one assigns a 
projection P u in the Hilbert space Ti such that TrjpP^} is 
the probability of measuring 10. This motivates the following 
generalization of quantum observables. 

Definition 2.1: [6, Def 3.1.1] Let D, be a set, & a er-field 
of subsets of f2, and Ti a Hilbert space. Then a "H-valued 
Positive Operator Valued Measure (POVM) on is a map 
A : & -> C + (H); E -> A{E) = A E such that 

1) ME) > M®) = for a11 E e,^; 

2) For any countable, mutually disjoint collection 
{E n } C T we have 



\ n / n 



where the series convergence on the right is in the 
weak operator topology; 

3) A(Vl) = Ifi, the identity operator on Ti. 
A POVM that corresponds to a physical experiment has 
a simple interpretation. The set ft is the sample space 
corresponding to experimental outcomes so that the cr-field 
& consists of the set of all events. The POVM A and a state 
p on Ti induce a measure p p ^(-) = Tr|pA(-)| on Q so 
that p p a(E) gives the probability of event E £ & . We use 
the notation A € E to denote the event that the measurement 
of POVM A resulted in a value in E £ 

Note that here il can have any general structure. This 
allows one to describe measurement apparatuses with out- 
comes that are physically e.g. multi-dimensional or on a 
manifold topology like the circle or sphere, which is not 
possible with standard von Neumann measurements. The 
latter are indeed equivalent to a special case of POVMs 
called Projection Valued Measures (PVMs), which require 
that f2 is a closed subset of the real line and that the 
range of A only consists of commuting projections. The 
unique correspondence between PVMs (A) and self-adjoint 
operators describing von Neumann measurements (A) is 
obtained through (1) by setting P^ = A(u>) for each X £ fl. 

It has been shown that any POVM on Ti can be viewed 
as the restriction to H of a PVM on a larger Hilbert space. 

Theorem 2.1: [6, Th.9.3.2] Let 51 be a compact metrizable 
space with Borel field & ', and A a POVM taking values in 
£+(%). Then there exists a Hilbert space K, D H and a 
PVM A K : & -> C+(JC) such that if P is the orthogonal 
projection from tC onto T-L then A(E) is the restriction of 
PA K (E)P to T-L for all Borel sets E. 
In principle, the existing filtering theory for PVMs [11], [12] 
thus covers the needs of POVM-based filtering, modulo a 
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proper lift of the Hilbert space. However, the latter is not 
unique and for reasons explained in Section I, it makes 
sense to look for a POVM theory that does not build on 
its reduction from a PVM. 

To generalize the second measurement postulate, 
Davies [6] introduces the notion of an instrument as a 
complement to POVMs. In the following we denote by 
CP(H) the set of all completely positive (CP) maps 3 
A : p n- A(p) on states p G S(H). 

Definition 2.2: [6, Def 4.1.1] Let O be a set, & a cr-field 
of subsets of S7, and H. a Hilbert space. Then a "H-valued 
instrument on fl is a map A : & — >• CP('H); E -> A E (-) 
such that 

1) A E > A = for all E G &\ 

2) For any countable, mutually disjoint collection 
{E n } C J- we have 

a (u„b„) = A En 

n 

where the series convergence on the right is in the 
weak operator topology; 

3) Tr{An(p)} = Tr{p}, for all p G S{H). 

If an experiment is set up so that the outcomes take values in 
some set f2, then for a quantum system initially prepared in 
state p G S(W) the measurement postulates for an instrument 
write: 

1) an outcome in the set £ C Q is obtained with 
probability P(E) = Tr lA^O) j; 

2) the state of the system conditioned on a measurement 
outcome in set E is A E (p)/V(E). 

Theorem 2.2: [6, Th.3.1.3, Th.9.2.3] If A is an instrument 
then for all E G & ', there exits a (non-unique) countable set 

{A n (E)} neN c C(H) such that 

A e(p) = ]T A n (E)pA n (E)* . (2) 

neiV 

Moreover, there exists a unique POVM A on fl associated 
to A such that for all E C £1 and p G S{T~L) we have: 

Tr{A B (p)} =Tr{pi £ } . (3) 

In particular, this unique POVM is given by 

A(E) = A n (E)* A n (E) for all E G & . 

Theorem 2.1 implies that it is also always possible to 
construct an instrument corresponding to a given POVM. 
However, there is no unique nor canonical way to choose 
the instrument without further information about the physical 
system. 

3 CP maps are the most general possible quantum evolutions, see [6, 
Sec. 9. 2] and [17, p. 251] for a discussion. The definition given below for 
instruments is different from that given in [6], wherein the assumption 
of complete positivity is replaced by positivity; see the discussion in [6, 
Sec.9.2]. 



A. System identification of POVMs 

Consider an experimental setup corresponding to unknown 
instrument A, associated to POVM A. In order to experi- 
mentally determine A, one can initialize the system being 
measured in some state p = |0) (0| and for any E G the 
probability of measurement outcome in set E is 

¥[AeE] = ^\A{E)\4>). 

With sufficiently many experimental outcomes, one can 
estimate the classical probability distribution (<j)\ A(E) \<f>) 
over all E C £1; doing this for different \<f>) and using 
polarization then allows to calculate A(E) itself. In order to 
ascertain the instrument A however, we must have access 
to the state A E (p) e.g. performing a state tomography 
experiment. This is often impractical in experimental setups. 

In the following sections, we use instruments in theoretical 
developments and to examine general properties of the 
measurement settings we consider. The goal is however to 
show that our final filtering equation works with POVM data 
only. 

B. Notation 

In the remainder of the paper we will use A, B, ... 
to denote instruments and A,B,... to denote the POVMs 
corresponding to these instruments. Also, if an instrument 
A corresponds to a von Neumann measurement, then we 
denote by A the associated self-adjoint operator. We will use 
the terms PVM and self-adjoint operator interchangeably. 

III. Commuting Instruments 

There is a very general notion of quantum conditional 
expectation that was first studied in [18] (cf. review in [19]). 
However, a definition of conditional expectation that is 
motivated by quantum measurements and is closely related 
to the classical Kolgomorov definition for random variables 
is discussed in [11], [12]. The classical definition of con- 
ditional expectation builds on joint probabilities, which are 
not obvious in the quantum context. Therefore the central 
idea in [11] is that in order to define a conditional expec- 
tation of two self-adjoint operators, the two operators must 
commute with each other. Indeed, in other situations the von 
Neumann measurement postulates imply that no consistent 
'joint' measurement can be defined. 

We now wish to generalize the [11], [12] definition of 
conditional expectation for two POVMs and for this we need 
to understand when it is possible to measure two POVMs 
simultaneously. In fact this depends on the commutativity of 
the instruments used to implement the POVMs. 

Definition 3.1: Suppose and (^2,^2) are two 

measure spaces and W is some Hilbert space. Then two H- 
valued instruments, Ai : -> CP("H) and A 2 : ,^2 -> 
CP("H) are said to commute with each other if for all 
E\ G &\ and E2 G ^2, there exist sequences of operators 
{A x n {E x ) : n G TVi} and {A 2 m {E 2 ) : m G N 2 } in C(H) such 
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that: 

• the instruments write (cf. Theorem 2.2) 



(4) 



){p) 



E /2 

nGJV 2 



A;;(£ 2 )p<(£ 2 )* for all p; (5) 



• for all m 6 AT^ and n € TVb we have: 

[>C(^i).^(^a)] - [^(^i).^K^»)1 = «• (6) 
Remark 3.1: For the special case of PVMs, Ni = N 2 = 1 

and the commutativity condition of Definition 3.1 is clearly 

equivalent to the commutativity of the associated self-adjoint 

operators A\ and A%. Note that in general, and even when 

satisfying Definition 3.1, a POVM instrument does not 

commute with itself (cf. proof of Theorem 4.3.1 in [6]). 

Now we consider the composition of instruments — the 
filtering application will involve one (actual) instrument 
on the probe and one (hypothetical, expressing our goal- 
variable) on the target system. 

Theorem 3.1: [6, Th.3.4.2] Suppose Aj, i = l,2,...n 
are instruments on some compact metrizable fij with Borel 
field J^. Then there exists a unique "joint" instrument A on 
ill X O2 X ... X il n such that for all Ei 6 ^ and p we have: 



L £iX_B 2 



?„(P) 



,(K.) 



, A 2 ,(E 2 )oA li(fil) (p). 



We now prove the first result of this paper. 

Theorem 3.2: Suppose Aj, i = 1, 2, . . . n are instruments 
on some compact metrizable f2j with Borel field and 
Ai are the corresponding POVMs. If the Aj are mutually 
commutative, then the POVM A corresponding to the joint 
instrument A is uniquely determined by the POVMs Ai , 
according to 

A(Ei xE 2 x ... x Ep) = Ai(E 1 )A 2 (E 2 )...A p (E p ). 
Moreover, the POVMs are mutually commutative, that is 

[Ai(Ei),M E 3)] = for a11 Ei> E h^ 3- 

Proof: We first prove the result for p = 2. From 

Definition 3.1 construct sequences of operators {A\ : n € 

N A } and {A 2 n : n e N B } satisfying (4),(5),(6). Select some 

events Ei € i^i, E 2 6 ,^ 2 and let E = E x x E 2 . Then for 

all p E S(H), we have 



Tr 



{pA E } = Ti[k E {p)) 
(Th.3.1) = Tr|A 2i(E3) o A 1)(Bl) (p)} 



(Th.2.2) 



Tr 



E 

meN 1 ,neN 2 



A n A mP A m A r 



ilr.ice pn.pem) Tr ^ p E A m A n A n A m 

meN 1 ,neN 2 



(Def.3.1) 



— Tr < p A 2 n A 2 n > 

{^1^2} , 



= Tr 



(7) 



where for notational convenience, we have written A l n for 
A l n (Ei), with i = 1,2. Taking the adjoint of the trace 



argument in the last expression and cycling p to the left 
finally yields Tr {pi^} = Tr |pi 2 ^i}- 

When adding a third instrument A 3 to this setting, we 
assume that A3 commutes with Ai and A 2 , but a priori 
the representations (2) used for these various commutations 
might be incompatible. We now show that the developments 
above can nevertheless be repeated. From Theorem 3.1 we 
know that A is uniquely determined by Ai and A 2 , and by 
Theorem 2.2 the POVM, A is uniquely determined by A. 
Therefore, given any two sequences of operators \B\ : n € 
Mi} and {B 2 : n e M 2 } such that 



K(EMP) = E B nP B n 

neMi 



(8) 



for i = 1,2 — here, we might choose {B^} that do not 
commute with the {B 2 } — we have 

E A n A%n = E B%n B ™ 



and also 
Tr 



{a e ( P )} = Tr{A l!(El 



) A 2,(£ 2 ) 



(P)} 



(definition) 



(put (8) in (7)) 



\ E b I b Ip b T b i 

Trip E B n B l E B ^ B l\ 
\ n£M 2 mSMi J 



This allows to prove the theorem for p > 2 with a recursive 
argument. ■ 

IV. Conditional Expectation 

We are now in a position to define the conditional expecta- 
tion of POVMs associated with two commuting instruments. 
Recall that we wish to find an expression for the conditional 
expectation that is expressed only in terms of the POVMs 
and not in terms of the instruments themselves. 

A. Basic definitions 

In this Section, we assume that % is a Hilbert space, Ha 
and 51s are two compact metrizable sets with Borel algebras 
J^A and ,^b, respectively, and A and B are two "H-valued 
POVMs on Ha and Hb corresponding to the instruments A 
and B. 

Suppose p is a state on H and / is any measurable function 
on Ha- Then the integral of / with respect to A over a set 
E E &a is defined by: 4 

f(u>)Tr{pdA(u>)}± [ f(cj)da A (u;) 

where the measure P p a{ ) — Tr|pA(-)| is a probability 
measure on (Ha,^a)- When a vector space is associated 

4 It should also be possible to define this integral as a Stiltjes integral over 
the measurable space £_)_(%) with the measure induced by the POVM. 
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to the POVM measurement results Ha, we can compute the 
expectation value of a POVM A in state p by: 



E p [A] 



jTr{pdA(u))\ 



Defining the inverse map A 1 : £ + (H) — > Ha of A in the 
standard way, we can rewrite this as 



E p [A] = / A~ (X)Tr {p dX} . (9) 

Now suppose A and B are such that A(,^ A ) C B{^ B ). 
Then for any E g ^4 we define 



^(w)) Tr (pdA(oj)\ 
(B- l A{uj)dn a (lj) 

4 p 1 ™ H <I '«' U p,i( jP «) 



if the limit exists. Here {P n } is some partition of Ha and Lu n 
is in the smallest set P n E such that -B(P„) > A(P n ). 

B. Main result 

For a given state p and two commuting instruments A, 
B with Ha a convex subset of a vector space, in analogy 
with the classical definition we would like to define the 
conditional expectation E[yl|i3] : Ha x — > of A 
with respect to B such that E[A|_B] is JFq -measurable an( j : 



u) A Tr{pd(AB)(Lj)\ 

A\B]((j)Tr{pd(AB)(u)\ (10) 



Ef 

for all E E ^ B . Moreover, referring to (9), we want to 
define a POVM E A \ B : H A -> associated to this 

conditional expectation, such that 



aMTr{pd(iB)(u;)} 
(e^.B(lj b )) Tr{pdB(u B )} (11) 



for all i? E The following proves that this is possible. 

Theorem 4.1: Suppose p is a state on % and A, B are 
commuting instruments with Ha a convex subset of a vector 
space. Then there exists a er-field Sf^ on Ha and a (p-unique) 
POVM E A]B : <S A -> £+{H) satisfying (11). We call E A]B 

the conditional POVM for A with respect to B. 

Proof: From Theorems 3.1 and 3.2 we know that the 
POVM AB on H = H A x H B is given by AB(E A x 
E B ) = A(Ea)B(E b ). Denote by p — p p Ab {-) me measure 
induced by AB on H. 

We now define two classical random variables 

a : (uj a ,u b ) 1 y io A and /3 : (uj a ,^b) i-> w b 

on O. Then we know from the classical Kolgomorov theory 
that the conditional expectation 

E[a|/3] : H -> A 



exists, is a (p-a.s.) unique random variable that is measurable 
with respect to the cr-algebra &p generated by j3 and satisfies 
for all E G J> 

a(cj)^(w) = / E[a|/3]d/x(w). (12) 

'w££ Ju£E 

Set Sfyj to be the cr-algebra generated by E[a|/3] in Ha- 
Clearly, &p = and for all E E&a the ^-measurability 
of E[a|/3] implies that E[a\j3]- 1 {E) £ & B . Therefore, we 
can define E^^ as follows: for any E E <Sa set 

-1 



E 



A\B 



(E)±B(E[a\0\ 1 (E) 



(13) 



Now (11) is simply a restatement of (12). Also, the p- 
uniqueness of E A ^ B follows from that of the classical con- 
ditional expectation E[a|/?]. ■ 

C. Example 

We now study a simple example to understand the ideas 
developed above. We note that it is possible to solve this 
example problem quite easily using some simple algebraic 
methods. However, the main utility of our novel formulation 
of conditional expectation becomes apparent in the case 
where the sample spaces of the POVMs are continuous 
valued as in the case of phase measurements [20]. 

In this example we assume that W s is a finite dimensional 
Hilbert space C d . We further assume that the system interacts 
with a single probe with state space H p and the system and 
probe are initially in the separable state p = p s ® p p in 
"Htot = H s (8)"Hp, where the initial state of the probe is a pure 
state p p = \tjjp) (ip p \. The system and probe undergo a unitary 
evolution U, after which a Hp-valued POVM B on H B , 
corresponding to an instrument B, is measured on the probe. 
We assume that B is discrete, that is there exists a discrete set 
{wi, oj 2 , ■ ■■} cH B such that B({w 1; w 2 , •■•})= Let A 
be a discrete H s -valued PVM on Ha = K, associated to the 
self-adjoint operator A on H s . Note that we have to associate 
Ha to a vector space, but we make no such assumption 
on H B i.e. for the measurement carried out on the probe. 
This will allow us later to consider probe measurements with 
results on manifolds, like the circle for phase measurements. 

For any W. tol valued POVM W we denote by W u the 
POVM U*WU on the same measure space as W obtained 
from the unitary evolution U on W (Heisenberg picture). We 
can extend any POVM defined on T~L S or T-L p to the tensor 
product "Htot = T~Ls ® H p ; making a slight abuse of notation, 
we use the same symbol for the POVM on the tensor product 
space. Clearly, any reasonable instruments associated to the 
^tot-valued POVMs A and B should commute, because 
those measurements are physically implemented on different 
tensor factors of H to i- Therefore, the instruments correspond- 
ing to the POVMS A u and B u also commute. Then the 
conditional POVM E Au ^ BU is well defined and the map 



u) B eE 



■E A 



v^Y^^b)) Tr{pdB u (ij B )} 



defines a real valued linear functional on the self-adjoint op- 
erators on H s that gives the expectation value of observable 
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A given measurement outcome B € E. Therefore, there 
exists a density operator p s on H s , which we identify with 
the post-measurement state of the system, such that 

Tr {A Ps } = <P(A). 

Using Equation (13) with the notation of the classical random 
variables a and f3 from the proof of Theorem 4.1, we have 
9(A) = J E K[a\j3]dfj,(cj), and for an elementary event 
E = {oJr} € &b some algebra gives us 

Tr {[p s ® p p ] B u (cj r ) (J2 QA OJ a A u (oj a )) } 

K{Ap s } = — ^ - ( ^ 

TT[[p s ®p P }BU( LUr )j 

Using the definition of (-) u , the cyclic property of the trace, 
and the fact that Tr {pA} = Tr j J2u a &Q A w a^( w a)} for a 
self-adjoint operator A, where CI a is the spectrum of A, we 
get 

Tr{u[ Ps ® Pp ]U* B(u r )A\ 

Tr{Ap s } = — L : 

Tr|f/[p s ®p p ]f/*BK)| 

Because this is true for all .4 we have 

Tt p {u\p s ®p p ]U* B{u r )\ 

Ps = — r ; r • 

Tr \U[p s ®p p ]U*B(u r )j 

V. Application to quantum filtering 

A. Measurement model 

Our measurement model is motivated from the discrete- 
time model used in [12] for filtering using PVMs. We con- 
sider a system with Hilbert space H s and a probe consisting 
of a sequence of subsystems n = I, 2, ... each with Hilbert 
space H n — H. So the probe is described on a Hilbert space 
H p = ®^—i% n , the combined state space of the probe 
and system is written Htot — Us ® H p . Suppose Cl n is a 
compact metrizable space 3 for n = 1, 2, . . ., let J£"„ a cr-field 
on fi n and B„ : & n — > CP(H n ) a H n valued instrument 
with corresponding POVM B n . Also, let (ft a, &a) be some 
measure space and A : CIa — > H s any system instrument 
with corresponding POVM A. We set H n i = H s for n = 
and recursively define H n ] = W n -i] ® for n > 1; 
similarly define f2„i = Ha x Oi x • • ■ x O n . Also set 

W(n = <8>i^„+l so that 'Htot = %n] ® W(„. 

Consider the initial state p' ot = \tpQ 0t ) (ipQ 0t \ on 'Htot, with 

10 = 1^0) ®~! IO • 

Thus \tp s ) is the initial system state and each \ip?) is the 
initial state of a probe subsytem on % n . 

We suppose that between time steps n and n + 1, the sys- 
tem interacts with the probe according to a unitary evolution 
operator U n on HgtHn, i.e. it interacts only with subsystem 
n of the probe. After this unitary evolution, the POVM B n is 
measured to be some u> n £ fi n . We wish to find a recursive 
equation for p n conditioned on the measurement outcome 
B n = u> n - 

5 lf f2 n is not compact then we can simply consider the 1 -point compact- 
ification of f2 n [6, p. 12]. 



B. Filtering equation 

With a slight abuse of notation, let U n : = Yli=i ^n> a 
unitary on H n ]. As in Section IV-C, instruments A and 
Bj commute for all i because they are defined on different 
Hilbert spaces in the tensor product space Htot, and clearly 
the same applies to instruments Bj,Bj for all i ^ j. Then 
the same commutations apply to the evolved instruments 

A(n) = U n] AU^ and B t (n) = U n] B t U* n] . 

Therefore, the conditional expectation 

E[A(n) \Bi (n)B2 (n) ••• B n (n)] is a well defined function 
of the measurement outcomes u>\, . . . ,u) n , and also the 
associated POVM 

%(n)|Bi(n)Bj(n)-i n (n) 

is well-defined. The main advantage of the commutativity 
condition is seen at this point: through the conditional expec- 
tation POVM, the conditional expecation of a system POVM 
A is an explicit function of the measurement outcomes if the 
instruments A, Bi, . . . , B„ are mutually commutative. 

Now we specialize to the case where the system in- 
strument A corresponds to a self-adjoint operator A : 
H s —> H s , i.e. the instrument A corresponds to a 
PVM.The linearity of the conditional expectation based on 
E X r m i> 1 \ i> 1 \ s / \ implies that 

A{n)\B 1 {n)B 2 (n)---B n (n) r 

* : / wTrl pl ot Ms, „A nfiM & , S[ 

/ 1 p A{n)\Bx(n)B 2 (n)— B n (n) J 

is a real valued linear map on the set of self-adjoint operators 
on H s - Therefore there exists a unique density operator p s n 
on H s that satisfies 

9(A) = Tt{p n A} 

for all self-adjoint A : H — > H. Therefore, p n is the esti- 
mated state of the system at time step n. Explicit expressions 
like (14) could be developed along the same lines as in 
Section IV-C. 

VI. Conclusion and Future work 

In this paper we show how quantum filtering can be per- 
formed in the POVM setting. We therefore formalize a notion 
of commuting instruments for two measurements, which 
gives a sufficient condition to define joint measurement 
statistics in terms of the associated POVMs only, without 
explicitly depending on the instruments. We then introduce 
the notion of conditional-expectation-POVM for two mea- 
surements with commuting instruments. On that basis, for a 
system-probe model, we establish the filtering equation for 
the system state conditioned on probe measurements with 
a general instrument, and we demonstrate that it is only a 
function of the probe POVMs. In a sequel to this paper [20], 
we will apply the developed ideas to derive a filtering equa- 
tion for discrete-time POVM phase measurements, which is a 
quintessential example of a POVM-type measurement where 
the associated instrument is not known. On a more general 
mode, we also intend to explore necessary conditions for the 
commutation of two instruments. 
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